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Abstract 

Let be a Wishart process of dimension 5, with values in the 
set of positive matrices of size m. We are interested in the large de- 
viations for a family of matrix-valued processes {5~ 1 X^\t < 1} as 5 
tends to infinity. The process is a solution of a stochastic differ- 
ential equation with a degenerate diffusion coefficient. Our approach 
is based upon the introduction of exponential martingales. We give 
some applications to large deviations for functionals of the Wishart 
processes, for example the set of eigenvalues. 
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1 Introduction 

Let B be a m x m matrix valued Brownian motion. We consider a Wishart 
process X t , solution of the following SDE, with values in <S+, the set ofmxm 
real symmetric non-negative matrices: 

(1.1) dX t = Jx t dB t + dB' t y/x~ t + 8I m dt, X = x, 

where x G and M' denotes the transpose of the matrix M. 
We recall the following existence theorem (see M.F. Bru pQ): 

if 5 > m + 1, and x G 5+ (the set of positive definite symmetric 
matrices), then (jl.lj) has a unique strong solution in 5+. 
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In fact, we can extend this result to a degenerate initial condition, and in 
the following, we shall allow x = 0. 

We shall look for a Large Deviation Principle for the 5+ valued diffusion 
with small diffusion coefficient: 



1.2) 



dXl = e(^/X] dB t + dB' t y/Xf) + 5I m dt, t<T 



X* = x 



with 5 > 0. For e small enough, according to the above existence result, (|1.2j) 
has a unique solution X\ e 5+ for t > 0. 

Note that this problem is equivalent to look for a LDP for the family of pro- 
cesses ( — X^ NS , t < 1) where xj: N ^ denotes a Wishart process of dimension 
N5, starting from Nx as iV — > oo. 

When m — 1, (jl.lj) is the equation for the squared Bessel process (BESQ) of 
dimension 5. 

In a previous paper [3], we studied large deviations for BESQ and squared 
Ornstein-Uhlenbeck processes. Note that the diffusion coefficient in the 
BESQ equation is not Lipschitz and the Freidlin-Wentzell theory doesn't 
apply directly (in the degenerate cases : x = or 5 = 0).We gave three ap- 
proaches; the first one was based upon exponential martingales, the second 
one uses the infinite divisibility of the law of BESQ processes (and thus a 
Cramer theorem) and the third method is a consequence of the continuity of 
the Ito map for the Bessel equation (not square), a property proved by Mc 
Kean [7j. 

We also refer to Feng [3] for the study of a LDP for squares of Ornstein- 
Uhlenbeck processes. 

In the matrix case, due to the restriction on the dimension 5, the laws Q 5 X of 
the Wishart processes are no more infinitely divisible. Moreover, we have no 
analogue of the Bessel equation for the square root of a Wishart process. 
Thus, we shall focus on the exponential martingale approach to extend the 
LDP in the matrix case. Since the delicate point is for a degenerate initial 
condition, we shall assume that x = 0. 

We denote by C o ([0, T]; <S+) the space of continuous paths ip t from [0, T] to 
<S+ such that ip = and ip t G for t > 0. 
The main result of the paper is: 

Theorem 1.1 The family P e of distributions of (X^;t G [0, T]), solution of 
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(|1.2j) . satisfies a LDP in C*o([0, T]; <S^) wift speed e 2 and good rate function 

(1.3) I i ( P) = \J ^(K^Ms)K(s))ds, <peC o ([0,T];S+) 

where k^s) is the unique symmetric matrix, solution of 

(1.4) k v (s)(p(s) + <p(s)k^(s) = 2((p(s) - 6I m ), s>0. 
Remark: In the real case (to = 1), we obtain (see [3]), 

8 Jo (p(s) 

The outline of the paper is the following. In Section 2, we prove an expo- 
nential tightness result for the distribution P e of X e . In section 3, we prove 
Theorem 1.1 using the approach of exponential martingales. In Section 4, 
we discuss the Cramer's approach, using the additivity of Wishart processes, 
when we put some restriction on the parameter 5. In section 5, we give some 
applications of the contraction principle to obtain a LDP for some functionals 
of the Wishart process. 

2 Exponential Tightness 

We follow the same lines as in |3J Section 2], that is, we prove exponential 
tightness in the space C a of a- Holder continuous functions with a < 1/2. 
Let a < 1/2 and set \\<p\\ a = sup 0<s _^ i<T ^r^pr where |.| is a norm on <S+. 
Since all the norms are equivalent, we shall choose a suitable norm and we 
consider in this section ||M|| = J2i<ij< m 

Proposition 2.1 The family of distributions P € of X e is exponentially tight 
in C a , in scale e 2 , i.e. for L > 0, there exists a compact set Kl in C a such 
that: 

(2.1) limsupe 2 lnP(X £ e K L ) < -L. 

£ >0 

Proof: Let us fix a' e (a, 1/2) and R > 0. The closed Holder ball B a ,(0, R) 
is a compact set of C a ([0, 1]). 
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Thus it's enough to estimate P(||X e || a ' > R). For simplicity, we assume 
T=l. 



\\X e \\ a i < \\M e \\ a , +5m 

where M e is the martingale defined by 

Ml = e(y^Q dB t + dB' t yfXt). 

Bounds for ||M e || a - We shall use Garsia-Rodemich-Rumsey's Lemma which 
asserts that if 



Jo Jo V p t-s / 



then 

•\t-8\ 



Ml - Ml\\ < 8 / V-\£K/u*)dp{u). 
Jo 



Take ty(x) = e ce 2x — 1 for some < c < 1/2 and p(x) = x 1 ^ 2 . So \& 1 (y) 
— log(l + y). This yields (see the same computations in [3]): 
(2.2) 



1 r l 



JO 



P [\\M e \\ a , > R < P / / exp C6-'V — wl ^)dsdt>K + l 



_ 2U Ml-Mt 



with K = \ [e^- K2 )- A - l) and K 2 = 2 sup u6[0jl] u 1 ' 2 ^' log ±. 
Now by Markov's inequality, 



(2.3) ^(imU^j^f/^hK 



It -si 1 / 2 



dsdt. 



Now, for a matrix M, 

exp(A||M||) = l[exp(X\M tJ \) < JJ[exp(AM y ) + exp(-AMy)] 

< m 2 max[exp(AMy) + exp(— A My) 

< m 2 ^ [exp (AM„- ) + exp ( - A My ) ] 

Thus, 

£[exp(A||M t e -M;||)] 

< m 2 ^ (P[exp(A(M^(t) - M/^s))] + P[exp(-A(M^) - M^s 

< 2m 4 maxP[exp(2A 2 (MM*)] 

ij 
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where we use in the last inequality the exponential inequality for continuous 
martingales 

E[exp(XZ t )} < E[exp(2X 2 {Z) t )}. 

Now, 



j/s 



e 2 / (XUu) + XUu))du 



< e 2 



Tr(X-)du. 



Set Y* := Tr(JQ), then, Y* is a squared Bessel process, solution of the 
following SDE 



(2.4) 



dY* = 2e v /T^ci/3„ + 5m (it 



with j3 a real Brownian motion. Thus, we obtain: 



E 



exp ce 



|m?-m;| 

|t-s|V2 



(2.5) 



< 2m 4 I E 



< 2m 4 



exp 



2c z e 



2 r -i rt 



(t-8) 



Yfdu 



1/2 



t - S 



E [exp (2cV 2 Y; £ )] du 



1/2 



(by Jensen's inequality). Thus, we obtain: 



2m 4 



1/2 



(2.6) P(\\M% >R)< — — \ sup E [exp(2c 2 e- 2 FJ 

K + 1 [«e[o,i] 

where X + 1 = C exp(oRe _2 /8) and C a constant. 
Now, 

E[exp(2c 2 6- 2 FJ)] = gr e " 2 [exp(2c 2 Xj] 

where Q£ denotes the distribution of a squared Bessel process, starting from 
x, of dimension p. The Laplace transform of the BESQ is known (jUj) and 
we obtain: for c < 1/2, 

Q^ e " [exp(2c 2 X„)] = (l-4cV) 2 . 

implying 

P(||M e || a , > R) < C m A m5e ~ 2 e- cRt ~ 2 / 8 
for a positive constant A. Thus, 

lim limsupe 2 lnP(||M e || a , > R) = -oo. □ 
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3 Proof of Theorem 11.11 

From the previous section, we need to prove a weak LDP, that is to prove 
the upper bound for compact sets. We assume that T = 1. According to [2], 
we shall prove: 

i) Weak upper bound: 

(3.1) limlimsupe 2 lnP(X e G B r (<p)) < -I(ip) 

where B r (ip) denotes the open ball with center G Cq ([0, 1];<S+) and 
radius r. 

ii) Lower bound : for all open set O C Cq ([0,1],; <S+), 

(3.2) lim inf e 2 In P(X e G O) > - inf I(cp) . 

e—>0 ipGO 

3.1 The upper bound 

We denote by Jvt m , resp. S m the space ofmxm matrices, resp. symmetric 
matrices, endowed with the scalar product: 

(A,B) = Tt(AB'). 

The corresponding norm is denoted by \\A\\ 2 . Set H = {h G C([0, 1]; S m ) : 
h G L 2 ([0, 1]; S m )}. For h G H let 

Mf = exp (±{J*Tr(h(s)(dX e a -5I m ds)) - l -(Z\ Z<) ,t < 1 

Z? = [ t Tr(h(s)^XjdB s + h(s)dB' s ^Xj). 
Jo 

{Z\Z £ ) t = 4 f Tr(/i(s)X s e h(s))cis. 

M e,/i is a positive, local martingale. In fact, using a Novikov's type criterion 
(see |9, Exercise VIII. 1.40], [3J), we can prove that M e ' h is a martingale, then, 



where 
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E{Ml' h ) = 1. 

By an integration by parts, we can write: 



e 2 



with 



and 



$(^5 h) = G(y?; h) - 2 [ Tr(h(s)ip(s)h(s)) ds 

Jo 

G(<p]h)=Tr(h 1 (<p 1 -8I m ))- f Tr((v9 s - 5sl m )h s )ds 

Jo 

for^eC ([0,l];5+). 

Remark: If (p is absolutely continuous, then, 

G(<p;h) = f Tr(h(s)^ s -5I m ds)). 
Jo 

For^GC ([0,l];<S+), /iGtf, 

(M e ' h * 

< exp(-^r inf $(^; h) ) E(M[' h ) 



w^hich yields : 



< exp ( -\ inf $0; /i) ) , 



limsupe'lnP(X e G 5 r (y)) < - inf 

e^O '4>eB r (tp) 

For h E H, the map y? — ► h) is continuous on Co([0, 1],«S+), so that 
limlimsupe 2 lnP(X e G B r (tp)) < -$(<p;h). 

e-,0 

Minimizing in h G H, we obtain: 

lim lim sup e 2 In P (X e G B r (tp)) < - sup $(</?; h). 

r^O e ^ heH 
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Proposition 3.1 For up e C o ([0, 1];5+), 

sup$(y?;/i) = I(y>) 

where I((p) is defined by ([1.3)1 . 

Proof: Since 93 G C ([0, 1]; <S+), £ Tr(h s (p s h s )ds > for ft ^ 0. Replacing 
by A/i for A G 1, we can see that 

J (p) := sup ft) = - sup - r - 



/ie// 8 heH j Q Ti(h s ip s h s )ds 

We assume that J (ip) < 00. We denote by ||ft||L 2 (v?) the Hilbert norm on 
C ([0, l];«5 m ) given by 



INIl 2 ( v ) = / Tr(h s ip s h s )ds. 
Jo 

The linear form G v : ft — > ^(99; /i) can be extended to the space L 2 ((p) and 
by Riesz theorem, there exists a function fc^, e L 2 ((p) such that 

(3.3) G(y?; h) = Tr(ft, a ^ s fc ¥ ,(s))ds 

Jo 

Thus, (y? is absolutely continuous and we have 

(3.4) / Tr(/i s (0 s - 5I m ))ds = / Tr(/i s <p s fcp(s))ds 



for all symmetric matrix h(s). Let k^ be given by ([1-4)1 . We refer to the 
Appendix for the existence of an unique solution of ([1.4)1 . Then, it is easy to 
see that ()3.4|) is satisfied for all h symmetric. Moreover, by Cauchy-Schwarz 
inequality, 

1 ,1 
TiihsUDsk^s^ds < I Ti(h s uj s h s ) 1/2 TT(k ip (s)uj s k ip (s)) 1/2 ds 

Jo 

< ([ Tr(/i^A)^) 1/2 ( [ Tr^V^s))^) 172 



with equality for h = k v . 

1 G 2 ( V ;h) 
Thus, - sup -— 2 = I{<p). 

Now, the equality between /(<£>) and J(y?) follows by density of H in L 2 ((p). 
□ 
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3.2 The lower bound 

In order to prove the lower bound, we first prove 

liminf e 2 \nP(X e e B r (<p)) > -I(<p) 



for all r > and for if in a subclass Ti of Co([0, 1]; <S+). Then, we shall show 
that this subclass is rich enough. 

Set TC the set of functions <p such that I(ip) < oo and s.t. k^ defined by (jl.4|) 
belongs to H . 

For <f G H, set h v = \k^. As in the previous subsection, we introduce the 
new probability measure 

p ■= M[' K P 

where P is the Wiener measure on C([0, 1]; A^ m , m ). Under P, 

B t = B t + ~ f (y/Tl h^s)) ds 
e Jo 

where B is a Brownian matrix on P. 
Thus, under P, X e solves the SDE 

dXl = e(y/X]dB t + d&ty/Xl) + {2{Xlh^{t) + K{t)Xl) + SI m )dt. 

Under P, Xf solution of 

d% = [2(%h v (t) + K(t)S> t ) + SI m ]dt 

i.e. 



1 

Since is continuous, this equation has ip as a unique solution; thus 



*t - «J/ m = 2{%h v {t) + M*)*t) = -(%k v (t) + k v {t)%). 



Xf — > (fit P a.s. 
and lim P(X e e B r (ip)) = 1 for every r > 0. Now, 

e > 

P{X*eB r {ip)) = P \X e E B r (<p)—^j- 



> exp | -i sup F(tP;\) ] P(X e E B r (<p)) 
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which yields : 



liminfe 2 lnP(X e € B T hp)) > - sup F(ip;\ 



and by continuity of F(.,h), 

lim liminf e 2 \nP(X e e BJip)) > Fi^hJ = I Up). 

r > e^O 

We now prove the: 

Proposition 3.2 For any p £ C ([0, 1]; S+) such that I hp) < oo, there 
exists a sequence <p n of elements of Tt such that ip n — >ip in C o ([0, 1]; <S+) 
and Ihp n ) — >Ihp). 

Proof: We follow the same lines as in the proof of the corresponding result 
for the scalar case in j^j. 

a) First, let us show that the condition I hp) < oo implies that 

lim ^ = 5I m . 

t — >o t 

From the scalar case, we know that: 

(3.5) lim = 6m , 

Indeed, Tr(Xf) satisfies a LDP (see (j2.4|l ) with rate function given by J(g) = 

I Jo ^gT^ ^ s an d J{g) < oo implies that lim t >0 ^ = ( see 0, jl])- 

From the upper bound, the condition I hp) < oo implies that J(Tr(</?)) < oo 
and thus ()3.5|) is satisfied. 
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Let us denote = Tr(\A\) and ||A|| 2 = (Tr(| A?)) 1 ' 2 for a matrix A. 

\\tp t -5tl m \\i = || / ((p a -6I m )ds\\i 



JO 

1 '* 



2 



s ||1 







< 2^ / ll^«M s )M s + / IIM^vJi^) 



o 



< / llv^TIb ||\/^M S )M S 







(Tr(^)) 1 / 2 (Tr(^( S )^^( S ))) 1 / 2 rf S 

/ /-t \ 1/2 / ,t \ 1/2 

< fy Ti((p s )ds\ I J Tr ■(k <p (s)(p a k (p (s))ds) 



Thus, 



'y-W m ||i< (ll^f^ ds ) ' (J^(K(^sK(s))ds^j ' 



According to ()3.5|) . the first term in the RHS is bounded and the second 
tends to as t tends to since 1(f) < oo. 

b) As a second step, we approximate ip by a function ^ such that fc^, G 
L 2 ([0,l];<S m ). Set 

VvW = 5t/ m , ^ < r/2 

if, r (t) = 8r/2 I m + (t - r/2)a r , r/2 < t < r 
A{t) = <p(t), t > r 

where the matrix a r is chosen such that ip is continuous in r. Let the 
solution of (jl.4j) associated with ip. Since k^(s) = for s G [0, r/2], and 
if)(s) is invertible for s > 0, k^ G L 2 ([0, 1]; iS m ). Obviously, r/v — in 

C ([0,l];5+). 

It remains to prove the convergence of I(ip r ) to I(ip), or that 

Tr(k^(s)ip(s)k^(s))ds — ► 0. 

r/2 

11 



pr 

:(k^(s)i/j(s)k-,p(s))ds = / Tr (k^(s)(ip(s) - 8I m ))d& 

Jr/2 

Tr (k^(s)(a r - 5I m ))ds. 



Jr/2 



Note that a r and k^ r (s) for s G [r/2,r] are diagonalisable in the same ba- 

sis with respective eigenvalues: (a^)i and kAs) = — - — fpr and that, 

according to step a), lim a\ r ^ = 5. Thus, for r small enough, 

r ► 

» r rr , (r) r\ 2 m 

/ Tr(k4sMs)k4s))ds = T^^^EK-^ 

A/2 ir/2 V <W 2 + (« - r / 2 K ^ 

and the last quantity tends to as r tends to 0. 

c) Thanks to b), we must find an approximating sequence p^ of p in TC 
for p satisfying k^ E L 2 . 

Let k^ be a sequence of smooth functions with values in S m such that k^> 
converges to k v in L 2 ([0, l],«S m ). Let p^> be the unique solution of 



( n ) 

V?o = 



Since 

ii Wii ^ 

\m II < 



f \\^ ) \\ds<2 / |bi n) || \\k^\\ds + 5, 
o Jo 

the Gronwall inequality shows that 



ii Hi i ^ 
sup sup \\ipl '\ \ < oo 

n te[o,i] 

where we have chosen the operator norm on the set of matrices in the previous 
inequality. Another application of Gronwall's inequality entails that: 

I I ( n ) I I n 

sup \\(pt-<Pt '\ \ — > 0. 



te[o,i] 



n— «x> 



Now, the convergence of I(p^) to /(<£>) follows from the convergence in L 2 
of /c^ to k v and the convergence in L°°([0, 1]) of p^ to p. □ 
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4 The Cramer theorem 

Let Q s x denote the distribution on (7(R, of the Wishart process of di- 
mension 8 > m + 1, starting from x G <S+. We recall the following additivity 
property (see 0): 

Qx ® Qy = Qx+y ■ 

Let 8 > m + 1 and take e = then X e , solution of (jl.2j) , is distributed as 
z27=i-^-i where Xj are independent copies of Q^. From Cramer's theorem 
(0 chap. 6), we obtain: 

Theorem 4.1 Lei 5 > m + 1. TTie family P e of distributions of {X^;t G 
[0, T]), solution of (1 1.2 J) , satisfies a LDP in Co([0, T]; wift speech e 2 and 
good rate function: 

(4.1) A*(<^) = sup (/ Tr(^) -A0i)V 

/ieA<([0,T],S m ) \</o / 



(4.2) A(/i) = In 



QifexpiJ Tr(X s dfi s )) 



The Laplace transform of the Qf, distribution can be computed explicitely 
in terms of Ricatti equation, extending to the matrix case, the well known 
result for the squared Bessel processes (see jH], [HI Chap. XI]). 

Lemma 4.2 Let \x be a positive -valued measure on [0, T]. Then, 

gj^exp(-| jT Tr(X a dj*.))) = exp(i Tr(F M (0)x)) exp(^ Tr(F M ( S ))d S ) 
where F^s) is the S m -valued, right continuous solution of the Riccati equation 



(4.4) F + F 2 = /i, F(T)=0. 

Proof: From Ito's formula, 

F^t)X t = F^(0)x+ [ F^(s)dX s + [ dF„{ 8 )X B 

Jo Jo 

= F M (0)x+ I F^{s)dX s + [ dfx(s)X s - [ F*{s)X s ds 
Jo Jo Jo 
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Consider the exponential local martingale 

Z t = exp Q J Tr(F IJt (s)dM s ) - 1 jT Tr (F^s)X,F^a))ds 
where M s = X s — 5/ m s. Then, 

Z t = exp Q ^Tr(F M (t)X t ) - Tr(F M (0)x) - 5 jf Tr(F„( S ))ds - jf Tr(X s d/i( S ))^ 

Now, X 4 is positive and F^(t) is negative 1 . Thus, TrpQi^i)) < and Z t is a 
bounded martingale. The Lemma follows from the equality E(Zq) = E(Zt)- 
□ 

Remarks: 

1. The condition F(T) = in (jOJ) is equivalent to F(T-) = -//({T}). 

2. Taking c?/x s = 20<5i(ds) where B is a symmetric positive matrix, we 
find that F^(t) = — 29(I m + 2(l -t)0) _1 , t < 1, from which we recover 
(seed): 

(4.5) 

Q s x (exp(- Tr^G))) = det(/ m + 26)^ 2 exp(- Tr(x(I m + 26)- 1 6)). 

For m = 1, this example is given in Subsection 8.3. 

Let us try to make the correspondence between ip and /i in (|4.1|) . If /i is a 
negative measure, then, from (|4.3|) . 

T /- T 1 5 f T 

TtfadiH) -AOi) = / Tr(<Mu t )--Tr(F_ 2/i (0)x)-- / Tr(F_ 2M (s))ci S . 

Since = — + F 4 2 ), an integration by parts gives: 

(4 T 6) 

jf Tr(<p t dnt)-AOj.) = lJ Tr(F.^(8)(<p 8 -6I n ))d8~ J Ti{F\{s)^ s )ds. 

The optimal function F(s) giving the supremum in ()4.6|) solves the equation: 

<f s - SI m = ip s F(s) + F(s)<p s 

that is F(s) = k (f (s)/2 where k v is the solution of (jl.4j) . and for this F, the 
RHS of fEED is exactly 
1 See the Appendix (A. 2) 
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5 Some applications 



From the contraction principle, we can obtain a LDP for some continuous 
functional of the Wishart process X e . 



5.1 The eigenvalues process 

Let (X e (t) = (Af(t), . . . , A^(i)); t G [0,T]) denote the process of eigenvalues 
of the process X e . 

Proposition 5.1 The process X e satisfies a LDP in Co([0, T], W?) with rate 
function: 



(5.1) j w = g£ 

Remark: (A e (t))t is solution of the SDE (see P 



from which we can guess the form of the rate function J in ()5.1|) since the 
drift b e in the above equation satisfies b e (X) — > 6. Nevertheless, since the drift 

b e (X) explodes on the hyperplanes {Aj = Xj} and the diffusion coefficient is 
degenerate, the classical results (see Theo. V.3.1]) do not apply. 

Proof of the Proposition I57T1 According to the contraction principle, 

J(x) = mf{I(<p); e.v.(ip) = x}. 

Write (ft = P^~ l K t Pt where A is the diagonal matrix of eigenvalues of cpt and 
Pt is an orthogonal matrix. Then, 

<p t = p t - l k t p t + p'r'AtPt + Pf'AtPt. 

We denote by k\ the matrix Ptk^fyPf 1 where k^ solves (II. 4j) . Then, 

Tr{k v (t)<p{t)k v {t)) = Tr{k t A{t)k t ) 

and 

kMXiit) + %(t)A,(t) = 2{Xi{t) - 5)5 l3 + R tj (t) 
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where the matrix R is defined by 



R{t) = PtPf'At + AtPtPf 1 - 
Now, it is easy to verify that Ru(t) = 0, thus: 

and the infimum of the above quantity is obtained for R = 0, corresponding 
to Pt independent of t. For this choice, I((p) = </(A) where A is the set of e.v. 
of if. □ 



5.2 A LDP for the r.v. Xf 

Proposition 5.2 The r.v. X{ satisfies a LDP, in scale e 2 , with rate func- 
tion: 

(5.2) K(M) = \ Tr(M) - i ln(det(M)) - ^ + ^ ln(5), M G <S+ 

Remark: For m = 1, 

K(a) = ~[(a-S)- ln(a/5)}, a > 

which corresponds (for 5 = 1) to the rate function obtained in the study of 
a LDP for a X2{n) distribution as n — ► oo. 

Sketch of proof: 

i) Since the application (p — ^(1) is continuous, we must minimize I(ip) 
under the constraint (p(l) = M. The optimal path (p solves the Euler La- 
grange equation (see [U], Chap. 7), given in terms of k v by: 

2^(s)+A;J(s) = 0,se(0,l). 

This leads to k~ l {t) = |J m + C and <p(t) = 5tl m + t 2 A with a matrix A 
determined by <p(l) = M. Note that this is the same path as in Section 4, 
Remark 2. 

Now, it is easy to verify that for ip(t) = 5tl m + t 2 (M - 5I m ), I ((f) = K(M) 
where K is given by (J5.2j) . 
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ii) Of course, we can compute K directly, using the Laplace transform 
(14. 5j) (with x = 0) and then, 

K(M) = sup{Tr(9M) + - ln(det(J m - 29))}. 
e 2 

The optimal 0o is given by M = S(I m — 2G>o) -1 . □ 

5.3 A LDP for the largest eigenvalue 

Let us denote by X e max the largest eigenvalue of the Wishart process X e . 

Proposition 5.3 The process {^ max (t),t 6 [0, T]} satisfies a LDP in Co([0, T 
with rate function given by 

Imax(f) = inf{J(a;),a; = (f,x 2 , . . .,x m ),Xi(t) < f(t) for i = 2, . . .m} 

where J is given by (j5.1|) . 

For f belonging to a class of functions T to be defined in the proof, 



(5.3) Imaxif) — 



-dt+ (m - 1) 







ft Jo f 



where f(t) = 5t + mf s < t (f(s) — 8s). 

Proof: We assume that the eigenvalues are given in decreasing order: Xi(t) > 
\ 2 {t) > ... > X m (t). 

According to the contraction principle, I max is given by the minimium of : 



m „t 



Xi{t) 



under the constraint {xi(t) < f(t),i = 2, . . . m} with X\ = f fixed 
Set 







y(t) 

F is a convex function on Co([0, T);R + ) and introduce the convex function 
G f (y)=y-feC([0,T);R). 

The problem is to minimize F(y) under the constraint G(y) < 0. 
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To /, we associate the measure fif associated to the Ricatti equation 
2(i f = H + H 2 on (0,T), H(T) = -2/i/(T) 

with H , = c/w-/). 

2/(t) 

Then, we define the measure djlf(t) = dfif(t) l(/(t)=/(t))- 

Let = {/; is a positive measure on [0, T]}. For / G JF, let us show 
that the Lagrangian 

L( k y,n) = F(y) + (G f (y),n) 

has a saddle point at (f,jlf), i.e., 

(5.4) L(l,fi)<L(l,fl f )<L(y,fl f ). 

for all y G Co([0, T);R + ) and all positive measure /i. 
The first inequality follows from 

(G f (l),ii)<0=(G f (l),jlf) 
since supp(fi f ) C {*,/(*) = f(t)}. 

For the second inequality, we must show that / minimize F(y) + (Gf(y),jlf). 
The optimal path of this problem of minimization solves the Euler- Lagrange 
equation (see [3]): 

< 5 ' 5 > l(iM = t ( " )+A/on(0 - n (t ( ">L = - Mn 

with g(a,b) = ^ 6 ~^ . The auxiliary function if t = ^y(7^ assoc i a ted to the 
optimal path y satisfies the Ricatti equation: 

2Jx f = H + H 2 ; H(T) = -2jl f {T). 

By the choice of /i/, it is easy to see that / solves the Euler-Lagrange equation 
()5.5p (or the associated Ricatti equation). 

According to Luenberger (Theo2, Section 8.4), the existence of this saddle 
point implies that : 

/ minimize F(y) under the constraint Gf(y) < 0. □ 

For a fixed time, we have the following result: 
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Proposition 5.4 The r.v. A^ iaa .(l) satisfies a LDP in M + with rate function 
given by 

& & & & 

(5.6) K max (a) = - - -ln(a) - - + -ln(5) if a > 5 



- - - ln(a) - - + - ln(5) j if a < 8 

The proof is immediate from (J5.2)) . We minimize K(M) under the constraint 
||M|| = a, where ||.|| denotes the operator norm. 



6 Appendix 

(A.l) On the equation AX+XA = B. 

Let A and B two symmetric matrices, with A strictly positive. We are looking 
for a symmetric matrix X, solution of the equation ( see (|l-4jl ): 

AX + XA = B (*) 

Since A is symmetric, let P and D be orthogonal and positive diagonal 
matrices such that A = P~ X DP. Then, from (*), the symmetric matrix 
X = PXP- 1 satisfies: 

DX + XD = PBP- 1 := B 

that is: 

B 

and Xij = — . Thus, X is uniquely determined. 

di + dj 

(A. 2) On the Riccati equation. 

We consider the Ricatti equation (see ()4.4|0 : 
(6.1) F + F 2 = /i, F(T) = 0. 
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or 

F(t) = C + /i(]0,t]) - / F 2 (s)ds 

Jo 

where C is chosen that F(T) = 0. We diagonalize F(t): F t = P^D t P t 
with D t the matrix of eigenvalues of F t and P t orthogonal. Then, the Ricatti 
equation can be written as: 

D{t) + D 2 {t) = p{t)mP-\t) + R t 

where R is a matrix, whose diagonal entries are zeroes. Set v = P/iP -1 , then 
v is a positive iS+-valued measure and the eigenvalues of F satisfy the scalar 
Riccati equation: 

di(t) + d 2 (t) = uu(t), di(T)=0 

with va a positive measure on [0,T]. We know (see jSj, Chapter XI) that 
di(t) < (di is related to the decreasing solution of the Sturm Liouville 
equation 0'/ = 0^). It follows that the matrix F(t) is symmetric negative. 
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